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Abstract 

We study the homogeneous but anisotropic cosmological models of Bianchi 
in presence of a massive scalar field <j> using the ADM Hamiltonian formalism. 
We begin to describe the main steps to find the ADM Hamiltonian of the General 
Relativity with a massive scalar field and then we study the dynamics of the flat 
Bianchi type / anisotropic Universe according to initial and final values of this 
Hamiltonian and sign of the potential. After a brief recall of the conditions neces- 
sary to isotropise an anisotropic Bianchi class A model with such a field, we extend 
them to a non minimally coupled scalar field by using a conformal transformation 
of the metric which casts the General Relativity with a scalar field into a scalar- 
tensor theory. The new line element then corresponds to the so-called Brans-Dicke 
frame, the former one being the Einstein frame. Then, we study the isotropisa- 
tion process of the Bianchi class A model when we consider the low energy form 
of the string theory without its antisymmetric tensor and the Brans-Dicke theory 
with some exponential or power laws of cj> for the potential. Finally, assuming an 
isotropic Universe such as all the metric functions behave as some power or ex- 
ponential laws of the proper time, we find the conditions such that the gravitation 
function and the potential of the scalar field are bounded as it is observed today, 
and compare them with the necessary conditions for isotropy. 

Please, review the English for style 

1 Introduction 

Scalar fields, which we will refer to by the greek letter <j>, are present in many parti- 
cles physics theories trying to describe the unification of the fundamental forces and to 
go beyond the interrogations of the standard model. Within this framework, they are in- 
terpreted like spin bosons. Hence the Higgs boson would make it possible to explain 
the mass of the particles while for supersymmetry theories, they represent additional 
degrees of freedom necessary to the existence of this symmetry. In this context, it thus 
appears logical to consider their existence at a cosmological level. Hence, they could 
be the source of the early times inflation or the recent acceleration of the Universe ex- 
pansion via what one calls a quintessent scalar field having a negative pressure. They 
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can also play a role in the variation of some constants of the nature like the gravitation 
constant or the cosmological constant. This one would then become variable and could 
thus explain the huge difference between its early times predicted value and its value 
possibly observed today. At a different scale, the scalar fields could also mimic the 
presence of a dark matter. Hence, in (DEL it was shown that they could explain the 
flatness of the rotations curves in the external regions of some galaxies. There are thus 
many reasons to consider the existence of some scalar fields in the Universe and it is 
what we will do in this work by studying the General Relativity with a scalar field or a 
scalar-tensor theory for which is non minimally coupled to the curvature. 

From a geometrical point of view, our present Universe is very well described by the 
homogeneous and isotropic models of Friedmann-Lemaitre-Robertson-Walker (FLRW 
models). Its expansion is then the same everywhere and in any direction. Such symme- 
try is absolutely extraordinary although it appears natural to us. To make a comparison, 
it is as extraordinary as a perfect straight line : mathematically it exists but in the real 
world no straight line is perfectly straight. There are several points of view to explain 
this. One can admit the isotropy and homogeneity of the Universe as an initial condi- 
tion. One can suppose the existence of a quantum principle which would select among 
all the possible geometries, the FLRW models or one can assume that initially our Uni- 
verse was less symmetrical than today and that it tends asymptotically to an isotropic 
and homogeneous state. It is this last point of view that we will adopt here by assuming 
that the Universe is initially homogeneous but anisotropic, which is the most immediate 
generalisation of the FLRW models. The homogeneous and anisotropic cosmological 
models were classified by Bianchi at the end of the nineteenth century in nine types 
listed by Roman letters going from / to IX. These types are divided into class A and 
class B, the Bianchi class A models being such as the trace of the structure constants 
is vanishing. 

In this work our goal will be to study the properties of the Bianchi class A models 
depending on those of the scalar fields with help of the ADM Hamiltonian formalism. 
In a first section, we give the main steps allowing to calculate the ADM Hamiltonian 
for the Bianchi class A models in General Relativity with a massive scalar field. In 
a second section, we show how the dynamics of the metric functions (contraction, 
expansion, etc) can be deduced from the initial and final values of the Hamiltonian 
and the sign of the potential for the flat Bianchi type / model. The isotropisation of 
the Bianchi class A models in presence of a minimally coupled and massive scalar 
field was already studied in papers |3]|4]. Hence, in a third section, using a conformal 
transformation of the metric, we show how to generalise these results to a scalar-tensor 
theory for which the field cf> is non minimally coupled to the curvature and study the 
cases of the Brans-Dicke theory and low energy string theory without its antisymmetric 
tensor when their potential are some power or exponential laws of <fi. Then, for these 
two theories, one remarks that when isotropy arises, the metric functions mainly tend 
to some power or exponential laws of the proper times. This leads us to examine what 
are the conditions, related to the parameters of these two types of laws, allowing the 
gravitation function and the potential of the scalar field to be bounded and to compare 
them with the necessary conditions for isotropisation. One discusses all these results in 
the last section. 
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2 The ADM Hamiltonian formalism 



There are three main Hamiltonian formulations of General Relativity : they are the 
Arnowitt, Deser and Misner (ADM) approach, the Dirac and the Kuchaf approaches. 
The ADM method consists in solving the primary constraints issued from the singu- 
lar Lagrangian of the theory and then developing the Hamiltonian formulation by using 
only the independent variables of the phase space. We have to remark that in the frame- 
work of the gauge theories, this method is far from being ideal : it indeed generally 
overlooks the covariance with respect to the Poincare group of symmetry and, in the 
case of constraints connected to a local invariance gauge, does not always succeed to 
highlight some aspects of the gauge invariance. In the present context, the ADM reso- 
lution of the constraints greatly simplifies the formalism (which does not include any 
more superfluous degrees of freedom) and the physical interpretation of the results (it is 
particularly interesting when one studies the spatial curvature effects on the singularity 
approach of the anisotropic models). 

The Dirac approach is issued from the Dirac constrained systems theory and in- 
cludes, without solving them, the constraints in the formalism ; it is particularly well 
adapted to the canonical quantification of the theory. Concerning the Kuchaf formula- 
tion, also interesting from the quantification point of view, it enlightens the geometric 
meaning of the General Relativity Hamiltonian formulation. In this work, we will fol- 
low the ADM method. 

The proof of the ADM results is particularly laborious ; these calculations are not 
often explicitly made in the literature. The B appendice of G. Rossi master thesis (For- 
malism hamiltonien en relativite generate et en cosmologie, Universite de Liege, Fac- 
ulte des sciences, Institut de mathematiques, 1973-1974), supervised by J. Demaret, 
shows the main technical steps. Some non published notes of R Tombal and A. Mous- 
siaux(Le formal isme hamiltonien en relativite generale (premiere version), Facultes 
universitaires Notre-Dame de la Paix, Namur) and a non published document of C. 
Scheen (Introduction au formalisme hamiltonien ADM de la relativite generale, Uni- 
versite de Liege, Institut d'astrophysique et de geophysique, 1992-1993), more system- 
atic and complete, show all the details of calculation. This section is inspired by these 
three works and proposes a summary of the main technical steps. 

The Hamiltonian formalism shows several advantages on the Lagrangian formal- 
ism. It allows writing the field equations as a first order equations system instead of the 
second order. The interpretation of the physical results is easier as shown by the clarity 
of the singularity chaotic approach of Misner using the ADM formalism with respect 
to the Belinskii-Khalatnikov-Lifshitz(BKL) one which uses the Lagrangian formalism. 

As a first step, we look for the Hamiltonian form of the General Relativity action : 



that we will generalise latter by taking into account a massive and minimally coupled 
scalar field. 

General Relativity is a typical example of a theory having the covariance properties 
for any coordinate change in spacetime : one says that it is parameterised a priori. In 
the classical Hamiltonian theory, it is possible to consider the time variable as a dy- 
namical variable ; the parameterisation of this theory reveals some constraints and the 
variational problem is to "extremalise" a form of the action where these constraints are 
introduced via Lagrange multipliers. Obviously, the constraints do not appear explicitly 
in the action Q. Moreover, in the Hamiltonian formalism, the time, separated from the 
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other variables, is considered as a parameter. The first thing to do is thus to rewrite the 
action Q by splitting space and time, i.e. by using a 3 + 1 spacetime decomposition. By 
doing this, we will see that the action can be written in the ADM Hamiltonian form : 



which will allow us to deduce the constraints. 

2.1 Writing the General Relativity action with the 3 + 1 spacetime 
decomposition 

2.1.1 3 + 1 spacetime decomposition 

The 3 + 1 spacetime decomposition consists of splitting it as a serie of spatial hyper- 
surfaces, parameterised by the t time. We start to define the lapse and shift functions. 

Be two hypersurfaces £(£) and £(t + dt) represented on the figure ^and their 
3-metrics, respectively ^gij(t, x k ) dx l dx 3 and ^gij(t + dt, x k ) dx % dxK Be the Pi 
point with coordinates (x % , t) on £(t). We define the P2 point as being the intersection 
of E(i + dt) with the normal to E(i) in Pi. The proper time interval dr = N dt 
between Pi and P2 then defines the lapse function N(xk,t). Let us define the P3 point 
on £(£ + dt) as being a point of this hypersurface having the same space coordinates 
as the Pi point. The P3 point coordinates are thus (x l ,t + dt) whereas these of P2 
are (x l — N z dt, t + dt). The vector binding P2 and P3 then defines the shift functions 
N t (xk, t). Be the P4 point on £(i + dt) with coordinates {x l + dx l , t + dt) and the P§ 
point on £(t) having the same space coordinates as P4, i.e. (x l + dx % , t). We define P5 
as being the intersection of the normal to E(t+dt) in P4 with E(i). The P5 coordinates 
are then (x l + dx l + N l dt, t). 

We are now able to express the line element ds 2 between the Pi and P4 points with 
help of the 3-metric ^gij, the shift and the lapse functions. Writing the Pythagoras 
theorem in the non Euclidean 4-space with signature (— , +, +, +), it comes : 




dt 




(4) g a/3 dx Q dx 

^ 9lJ (x k ,t)(x l (P 5 ) - x l (Pi))(xi(P 5 ) - x?(P x )) - dr 2 
(3) g tj (x k , t)(dx i + N' 1 dt){dx j + N j dt) - N 2 dt 2 



from which we get for the metric : 




or, putting AT, = ^ g t j N 3 : 




(3) 



and using (4) 5 Q /3 (4 V 7 = $Z '■ 




(4) 
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P4(xi+dxi,t+dt) 



P5(xi+dxi+Ni dt,t) 
• P6(xi+dxi,t) 
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Fig. 1 — The 3 + 1 spacetime decomposition. 



To calculate the determinant of the 4-metric, we use the Frobenius-Schur the 
orem which shows that if A, B, C and D are four square matrices, the matrix determi 
nant : 

'A B 



is det(A) = det(.D) det(A - BD^C). We then deduce : 



V-^g = N^)~g (6) 

Before continuing, we have to define the concept of extrinsic curvature. The ex- 
trinsic curvature characterises the way in which a variety is included in a space of 
higher dimension. For example, a paper sheet with two dimensions which one twists 
in a space with three dimensions has an extrinsic curvature with respect to this space. 
On the contrary, our Universe has an intrinsic curvature which does not need some 
higher dimension to be defined. In the case which interests us, the extrinsic curvature 
of a spatial hypersurface is a measure of the direction variation of the normal n to the 
hypersurface £(i) between two infinitely close points on E(£) and is defined by : 

Kij = -n i;j = -ej ■ ViH (7) 



2.1.2 The Gauss- Weingarten and Gauss-Codazzi relations : rewriting the action 

When the extrinsic curvature is known, one can express the covariant derivative 
of the basis vectors ej of the spatial hypersurface £ in the spacetime, depending on 
quantities only related to E. These are the Gauss-Weingarten relations which write : 

< 4 'V t e>^, 3 u+< 3 )^4 (8) 

The Gauss-Codazzi relations try to express the spacetime intrinsic curvature de- 
pending on the intrinsic and extrinsic curvatures of the hypersurface. They express as : 

Kik\j - K ij\k (9) 
-(KijK m k - K ik K mj ) + ^Rmijk (10) 
K iktn + K^K im (11) 

where | stands for the covariant derivative on the hypersurface. These relations allow 
rewriting the scalar curvature depending on the intrinsic and extrinsic curvatures on the 
hypersurface. Indeed, one can rewrite : 

WR = - 2^Ri j0 

The right hand side of this equation only contains some terms expressed by the 
Gauss-Codazzi relations and if we define Tr(K 2 ) = K^ k K ]k and K = K\ = Tr(K), 
one gets : 

( 4 )i?*j = g ik g^R kmij = - Tr(K 2 ) + K 2 

and : 

(4) <-o = 9 jH4) Rk0j0 = K, n + Tr(K 2 ) - K kj g% 
Since moreover, one can show that : 

g jk = 2R jk 



J hjk 
^ ^ Rmijk 
^RiOkO 
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it comes for the scalar curvature : 

(4) R= {3) R + Tt(K 2 ) + K 2 -2K n (12) 

Using this last relation and 0, we rewrite the Hilbert action Q in the following 

way : 

S= [ Ny/Wg~( {3) R + TT{K 2 )- K 2 )d A x-2N [ Ky^d 3 x (13) 

Jm JdM 

The surface term may be removed by imposing specific conditions on the variety 
border or by adding to the initial action another surface term compensating the one of 
(1131 - in this last case, we remove the surface term by taking as the initial action : 

S= [ ^R^^g~d i x + 2 f K^/^h,d 3 x (14) 

J M JdM 

In the case of a closed space, removing the surface term is not a problem : there is 
no consequence on the field equations when one varies the geometry inside the surface 
border of the variety. However, for open spaces which are asymptotically flat, it is 
necessary to add a surface term. 

Whatever the way in which one gets rid of the surface term, one obtains : 

S= [ N^/(Vg~(^R + Tr(K 2 )- K 2 )d 4 x (15) 
Jm 

In what follows, we are going to rewrite the action J 141 with help of the 3 + 1 
decomposition 

2.1.3 Writing the action as a 3 + 1 decomposition 

To write the above action as a 3 + 1 decomposition, we have to make this transfor- 
mation for : 

1 . the Christoffel symbols which are : 

r I/3 = 2~9 l5 {gas,p + gsp, a - g a p,s) 

2. the Ricci tensor which writes : 

(4) p pa pa , pa p/3 pa p/3 

y — ij,ot 1 ia,j ~r 1 ij L a/3 1 ja 

3. the scalar curvature which is ^R = ^'R°, and thus the action. 

To rewrite the Christoffel symbols V, we introduce the following definition : 

N j 
5 N 

and define the components of the tensor A, the Christoffel 3-symbols relative to the 
hypersurface, as : 

\? -^-(3)ajm((3) , (3) (3) \ 
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After tedious calculations, we get the following forms for the Christoffel symbols : 



1 00 






pi 

1 00 






L ik 






r o 

1 Oi 


Nu 

— Li _ k P 




1 i0 







We inject them in the Ricci tensor spatial components ^'Rij which express de- 
pending on the Christoffel symbols as : 

(4) n pa pa , pa p/3 pa p/3 

"+3 ~ 1 ij,a 1 ia,j 1 ij 1 a/3 1 i(3 L ja 

Then, it comes : 

+ A% A - + If^g + X^* + i k K ij]h (16) 

where ^Rij is the Ricci tensor of a hypersurface and is thus conventionally defined as 
a function of its Christoffel symbols as : 

(3) p \k \k I a k A ( A & A ' 

fi iJ — iv ij,fc — 2 Hk,j T ^ij^kl ~ £ Hl lv jk 

To continue our calculation, we should also rewrite explicitly the ^'Roi ar, d ^-Roa 
components of the Ricci tensor. However, these calculations are definitely harder than 
those which lead to the purely spatial components d!6l > of the Ricci tensor. In fact in 
General Relativity, it is enough to use a reference frame which simplifies the calculation 
without occulting the information related to the reference frame degrees of freedom. 
We will use the reference frame defined by the relations : 

_ 1 c) N l d 

n ~ Ndt ~ ~N~dx* 
_ ^ d 
1 dx 1 

In this particular frame, we have g nn = n-n = —1, g n { = n- e*j = 0, <7y = e*i • ej. 
The scalar curvature writes then as ^R = 2G™ + 2g^^Rij, where G a fj is the 
Einstein tensor. In virtue of the Gauss-Codazzi relations i ll Oi l, we get : 

G" = -i( (3) i? + K 2 - KijRv) 

We introduce this expression and (1161 in the scalar curvature ( 4 ) R and get : 



ZKijK** +4Kitf i + 2g«K ij[k £ k 



8 



where ^R is the scalar curvature relative to the hypersurface ; it expresses as : 



S 



Finally, the action in the 3 + 1 decomposition will take the following form : 

8K U dK 



/ d 4 x y/g 
Jm 



2JVjj + ZN^Kvtf + 2NK ij ^ i 



(17) 



2.2 Identification of the action 3 + 1 form with its form in the 
Hamiltonian formalism 

In this section, we wish to show the equivalence between the form dl7l of the action 
and its form adapted to the Hamiltonian formalism : 



S = I d A x 

M 



For that, we define the canonically conjugate momenta 7r u : 

* ij = y/g(g ij K-KV) 
and introduce the superhamiltonian as : 

Co = -( (3) R + k 2 - K^) = -Vg (3) R + M K n Rij - k 2 

However, we calculate that : 



KijK 13 - K 2 



fJ 



Tr(7r 2 )--(Tr(7r)) ; 



and consequently, it comes : 



c = -V9 {3) R + ^p 



Tt(7r 2 ) -(Tr(Tr))' 



In addition, we can also show that 

dir ij 



dt ~ V9 [ 9 at m 



Last, we introduce the supermomenta Ci via : 

= -2N\K\ - SiK^y/g = 2N* [-y/g{K> - SfK)] 
that we can rewrite : 



1.7 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



C i = -2ni {j = 2y/g(Ki-6iK) lj 
C l = -27rg = 2Jg(K* - g^K) b 

By inserting all these results in the action dl8> . we recover, after some additional 
calculation, the action written with the 3 + 1 decomposition. 
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2.3 Formulation of the ADM constraints in the General Relativity 

The canonical momenta tt ij are naturally defined by : 

^ ■ 6L 

TT = — 

L being the Lagrangian. The action of the Hamiltonian formalism writes : 

S = j (tt« ^-NC ~ N l C t ) d*x (24) 

The demonstration of the equivalence between the above expression and the form 
of the action is similar to assure the equivalence between the actions d!5l > and 0. 
Consequently, by varying i24l with respect to the lapse and shift functions, considered 
as some Lagrange multipliers, we get the constraints : 



G) = -Vs (3) i?+^ 



Tr(7T 2 )-i(TrW) 2 



(25) 

y/!J L z J 

d = -2tt^. = (26) 

The constraints rule the dynamics of the geometry and are, in the same time, some 
initial values conditions. Because of the constraints C M = 0, it is possible to choose 
freely the fields ^'g%j and 7r y on the initial hypersurface £(io)- The dynamical equa- 
tions rule the changes of the intrinsic geometry and the extrinsic curvature of a hyper- 
surface when one moves from a hypersurface to a neighbour hypersurface. 

If the constraints are satisfied on E(£o) and if the canonically conjugate fields 
evolve and check the dynamical equations, then the constraints are conserved all the 
time. 

In virtue of the constraints, four of the ( 3 '<7-y and 7T y fields may be expressed with 
the others ; moreover, imposing the coordinate conditions fixes four of the remaining 
fields. Only remain two couples of canonical variables ; the following paragraph shows 
how the ADM approach is brought back to these two physical degrees of freedom. 

2.4 ADM formulation of the General Relativity with a minimally 
coupled and massive scalar field 

To know how are modified the constraints equations (125b and i26\ when a scalar 
field <fi is present, we are going to write the 3 + 1 decomposition of the part of the action 
with the scalar field : 



S = 



V-^ffrf 4 x (27) 



where ui(4>) and U((j>) are respectively the Brans-Dicke coupling function describ- 
ing the coupling between the metric and the scalar field and the potential of the scalar 
field describing the self coupling of <p. First, we rewrite the term containing oj in the 
following form : 



-(3/2 + uj^^^V^g = (3/2 + w)^ 2 ^ 2 ^" 1 ^ (28) 

Taking into account this last expression, we can express the scalar field conjugate 
momentum as : 

dl 

tt = -r = (3 + 2uj)<j )( / ) - 2 N- 1 ^ (29) 
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where / is the Lagrangian of the above action. From this last equation, we deduce : 

N (f> 2 



^'^g3 + 2u 
which allows us to rewrite J28t in the following way : 



(30) 



-(3/2 + ^)<K<^»V (4) g= ~' t-tV9^ 



N 2 



2 /nn_ (3/2 + u) i 

" " ~ ^ iV vy "*T(3 + 2c)2 
1 4? N a 



+ ^ 

To the constraint (1251 . we have thus to add the terms Cq^ coming from the Brans- 
Dicke coupling function and the potential, Cq0 being such as : 



Thus, using the expression for it^ given by J29i and of <j), it comes : 
1 6 2 1 



(31) 



Co 



2 3 + 2oj ^fg 
1 6 2 1 



^77 + C/ V5 



(,2 



1 



2 3 + 2^^5 91 3 + 2w^<? 
1 6 2 1 



2 3 + 2lo Jg 

The final form of the Cq constraint, taking into account the scalar field, is thus 

1 6 2 1 



C = -^R+^- 



Tr(7r 2 )-i(TrW) 2 



2 3 + 2w y/g 



< + Uy/g 



The general form of the diagonal homogeneous but anisotropic metrics of Bianchi 
is defined by : 

ds 2 = -iV 2 ^ + i?2 e -20^/3 + + V^/3- (LL) l ) 2 +e /3 + -V^/3- (tJ 2 ) 2 +e -2^ +(tj3) 2^ ^ 

where the uji are the 1 -forms specifying each Bianchi model. This metric is anisotropic 
because it describes a Universe where the expansion rate depends on the direction. The 
Misner parameterisation J3]|6| of the conjugate momenta expresses as : 



Pk = 2?r7r fc - oir&kKl 



6pij = diag(p+ -I- V3p_,p+ - V3p-, -2p^ 
Then, we can rewrite the constraint on the following form : 



3 -3fi 



Mr- 



k\2 



-M) 



24tt 2 



(pI+p 2 -) 



2i?3 e - 3 ° (3 + 2u) 



Rle~ 3n U 
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7r being the number in this last expression, and the action becomes : 

S = J p+ df3 + + p_ d(3- + dcj) - H dfl 

with the ADM Hamiltonian H = 27T7t^ and the scalar field conjugate momentum 
Ptp = ttH^,. By using the constraint Co = 0, we then deduce for H : 

H 2 = p 2 + + v 2 - + 12 > + 24n 2 R 6 Q e- 6n U + V(tl, /?+,/?_) (33) 

where V(£l, f3+ , /3_ ) are the curvature terms specified in the table^for each Bianchi 
class A model. For the class B, the ADM Hamiltonian formulation does not give the 
correct field equations as explained in 1 7 1 and we have to redefine the divergence theo- 
rem in a non-coordinated basis to solve this problem and get a correct variational prin- 
ciple. The physical degrees of freedom have been isolated, but with this form the theory 
is no more covariant - the constraints have been solved and the coordinate conditions 
have been fixed. The loss of covariance is obvious besides : the ADM Hamiltonian is 
not null, while the cancellation of the Hamiltonian is characteristic of the constrained 
systems. 



Bianchi type 


V(Sl, /?+,/?_) 


I 





II 




VI , VII 


247r 2 i?4 e - 4n + 4 ^+(cosh4V3/?_ ± 1) 


VIII, IX 


24 7 r 2 i? 4 e- 4n [e 4 ' 3 +(cosh4V3/3^ - 1) + 
l/2e" 8 ^+ ± 2e" 2/3 + cosh2V3/3_] 



TAB . 1 — Form of the curvature potential V (O, 0+ , 0— ) for each Bianchi model. 



3 Dynamical study of the General Relativity with a mas- 
sive scalar field for the flat Bianchi type / model 

The results of this section are issued from papers 1 8 1 published in Classical and 
Quantum Gravity copyright 2000 IOP Publishing Ltd. We aim to describe the global 
properties of the metric functions (expansion, contraction and extrema) depending on 
Hamiltonian initial and final values and the sign of the potential when we consider 
the General Relativity with a massive scalar field in a flat Bianchi type / model, the 
corresponding Hamiltonian and metric being respectively defined by J33I and i32\ with 
(wi, cl>2, UJ3) — (dx, dy, dz). From [9] 1 , we know that the lapse function N takes the 
form : 



127ri?ge 



-30 



N = ^ (34) 

and then using the relation dt = —Ndil, it comes : 

dt ~ 2Ho[ dt d t )€ ~ 2H ° e HN 

'(see p!830 for a detailed calculus) 
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To obtain the global properties of the metric functions, we have to study the sign of 
the quantity J35i which depends on H — p^ . For sake of simplicity, we will consider 
a potential with a constant sign. We will see later how to extend our results when the 
sign of the potential varies. For the Bianchi type / model, the Hamiltonian equation 

where the dot means a derivative with respect to ft, shows that H and H are monotonic 
functions of il with constant sign. Then, from the lapse function, it comes that is a 
monotonic function of the proper time t and so for H (t). If t varies from tj„j to t en d, 
H will thus evolve monotonically from H(ti n i) = Hini to H(t en d) = H en d. 
Moreover, from the Hamiltonian equations 

i>± = -^ = (37) 

we derive that the conjugate momenta p^ are some constants (this is not the case in 
presence of curvature). Hence, the equation d35l may only have one zero and when the 
potential has a constant sign on a proper time interval, there is to the more one single 
extremum for the metric functions. 

From J35I . it is clear that the metric function derivative will vanish if the three con- 
ditions C\, C2 and C3 are true : 

- Cy : H and p^ have the same sign 

- C2 and C3 : p^ belongs to the interval defined by ffj n j and H en d 

Hence, we have to consider the following four cases for which we give the variation of 
the metric function with respect to the proper time t : 

Case la : U < 0, H and H > 

Taking into account J34i and the relation dt = —Ndft, it comes that when O increases, 
t decreases. The Hamiltonian is thus a decreasing function of t and the three conditions 
Ci write : 

- d : p i:j > 

- C2 : Hf in — pij<0 

- C3 : H ini — Pij>0 

Whatever the case we consider, if C2 or C3 are false, respectively C3 or C2 are true. In 
addition, in the present case, if C3 is false, C\ is true. 

If the three conditions are true, the metric function has a maximum in the proper 
time of the Einstein frame since the Hamiltonian will be equal to p^ for a value of O. 

If C\ is wrong, it is increasing since then H — ptj has always the sign of H. 

If C2 is wrong, the Hamiltonian is always larger than py, and the metric function 
is again increasing for the t time. 

If C3 is wrong, C\ is true, and the metric function decreases since the Hamiltonian 
is always smaller than pij . 

The same reasoning will hold for the other cases. 

Case lb : U < 0, H and H < 

When fl increases, t is increasing. The Hamiltonian is a negative and decreasing func- 
tion of these times coordinates. When C2 is wrong, C\ is true. The three conditions can 
be written as : 
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- d : pij < 

- C 2 : Hf in — pij<0 

- C3 : H ini — Pij>0 

If they are all true, the metric function has a maximum. 

If C\ or C3 is false, it is decreasing. 

If C\ is true, C2 is false and it is increasing. 

Case 2a : U > 0, H < and H >0 

When increases, t decreases. The Hamiltonian is a positive and decreasing function 
of f2 and then an increasing function of t. When C2 is wrong, C\ is true. The three 
conditions can be written as : 

- Ci : > 

- C 2 : Bend - Ptj>0 

- C3 : H ini — pij<0 

If they are all true, the metric function has a minimum. 
If C\ or C3 is false, it is increasing. 

If C2 is false, C\ is true, and the metric function is decreasing. 

Case 2b :U > 0, H > and H <0 

When £1 increases, t increases. The Hamiltonian is a negative and increasing function 
of the two time coordinates. When C3 is false, C\ is true. We obtain for the three 
conditions : 

- Ci : p^ < 

- C 2 : H end - p. y >0 

- C3 : H ini — Pij<0 

If the three conditions are true, the metric function has a minimum. 

If C\ or C2 is false, it is decreasing. 

If C\ is true, C3 is false, the metric function is increasing. 

All these results are summarised in table [2] They can not be extended to the other 
Bianchi models with curvature. They will be commented in the discussion. In the next 
section, we study the isotropisation of all the Bianchi class A models for a non mini- 
mally coupled and massive scalar field. 







H,H > 0, 


H, H, 


H, U > 


H, U > 0, 






U < 


U < 


H < 


H < 


Cx, 


C2, C3 : true 


Maximum 


Maximum 


Minimum 


Minimum 


Ci 


false 


Increasing 


Decreasing 


Increasing 


Decreasing 


c 2 


false 


Increasing 






Decreasing 


c 3 


false 




Decreasing 


Increasing 




c 2 


false, C\ : true 




Increasing 


Decreasing 




c 3 


false, C\ : true 


Decreasing 






Increasing 



Tab. 2 - Dynamical behaviour of a metric function in the proper time depending on 
the signs of the potential, the Hamiltonian and its initial and final values. 
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4 Bianchi class A isotropisation in scalar-tensor theory 



The isotropisation of the Bianchi class A models in General Relativity with a mas- 
sive scalar field was already studied in several papers 1 3 4|, using the ADM Hamilto- 
nian formalism, when the Universe expands forever (fl — > — oo). Defining the function 
I of the scalar field : 

f W± 
C/(3 + 2w) 1 /2 

these papers rested upon the following assumptions : 

1. The potential is positive. This is reasonable if we suppose that it could play the 
role of a varying cosmological constant. 

2. 3 + 2uj is positive, i.e. the scalar field weak energy condition is respected. 

3. Near the isotropic state, the perfect fluid density parameter f2 m , the scalar field 
energy density parameter Sl^ and the shear parameters £± tend to some equi- 
librium values, being non vanishing for f2^, This assumption physically defines 
what we called a class 1 isotropisation 2 and is in agreement with what we observe 
today. For a more precise mathematical definition, see for example 1 10 1. 

4. The isotropic state is reached sufficiently quickly. Physically it means that the 
variations of E±, fl m , fi^ and I in the neighbourhood of the isotropic state 
compared to their constant equilibrium values do not have any dynamical con- 
sequence on the asymptotical behaviour of the metric functions and potential. 
Among others, when in the neighbourhood of the isotropy I 2 tends to a constant 
£\ with a small variation S£ 2 , one assumes that J £ 2 + 5£ 2 d£l — > £q^I, the effects 
of S£ 2 thus being negligible. Anew, for a more precise mathematical definition, 
see itTPl . 

With these assumptions, we showed the following results concerning the isotropisation 
of the Bianchi class A models : 



Let us consider the General Relativity with a massive scalar field and the function 
I. The asymptotic behaviour of the scalar field near the isotropic state is given by the 
asymptotical solution of the following differential equation when Q — > — oo : 

A necessary condition for a class 1 isotropisation is that £ 2 tends to a constant such 
as £ 2 < 3 for aflat model, 3 being the spatial dimension, or £ 2 < 1 for a model with 
curvature, 1 being the spatial dimension minus 2. If £ tends to a non vanishing con- 
stant, the metric functions tend to t and the potential disappears like t~ 2 . If £ tends 
to vanish, the Universe tends to a De Sitter model with a cosmological constant. 

In what follows, we will use these results to study the isotropisation of the Bianchi 
class A models for a scalar-tensor theory, i.e. with a scalar field non minimally coupled 
to the curvature. The action of this theory can be deduced from the action (I27> by using 
a conformal transformation of the metric 

5q/3 = G~ 1 g a/ 3 



2 There are three isotropisation classes) 101 . The class 2 is similar to the class 1 but with — ► near the 
isotropisation state. For the class 3, the density parameters Qj, and/or Q m do not reach some equilibrium 



values and may oscillate 1 1 |. 
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where G is the gravitation function depending on the scalar field. The barred quantities 
are the ones of the new frame known as the Brans-Dicke frame, the former one being 
called the Einstein frame. Using this metric, the action \21\ is rewritten in the form : 




(39) 



with the following relations binding the quantities of the two frames : 



dt 



Ndn = VcF^dt 



(40) 



e 



(41) 



U = UG 



a 



(42) 



uj + 3/2 
02 



3/2(G- 1 )lG 2 +u>G<f>- 



-l 



(43) 



The necessary conditions for isotropisation of the Bianchi class A models with a non 
minimally and massive scalar field remain identical to those found in the Einstein 
frame. The conformal transformation does not change the fact that the Universe isotropises : 
if it occurs in the Einstein frame, it will occur in the Brans-Dicke frame. To express 
these conditions in the Brans-Dicke frame, it is enough to replace U and uj by then- 
expressions ( I42I43I > in the £ function. On the other hand the asymptotic behaviours of 
the metric functions and the potential can not be expressed any more in a closed form 
according to the proper time i because of the differential relation J40i which implies 
an integral of the gravitation function G. 

The two following subsections present two methods to study the Bianchi class A mod- 
els isotropisation for the scalar-tensor theory specified by the action d39l on the basis 
of the above results. 

4.1 Isotropisation of a scalar-tensor theory when the functions G, 
uj and U are known 

Supposing that one knows G((f>), U(<p) and oj(4>), we will use the following algo- 
rithm to calculate the asymptotic behaviour of the metric functions when the Universe 
isotropises in presence of a non minimally coupled scalar field : 

1. We calculate (f>(fl) using J38i what enables us to evaluate £ and thus the asymp- 
totic form of e according to the proper time t. 

2. We thus get n(t), cf)(t) and then t(t) thanks to d40l . 

3. Using d4ll . we obtain e~°. 

Let us apply this method to the Brans-Dicke theory and the string theory at low energy 
with power and exponential potentials of the scalar field : 

- Brans-Dicke theory with a potential in power of the scalar field 
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This theory is defined by G 1 = </5, u> = uiq and U — (j) k . Consequently U = <fi k 2 
and £ is a constant being 

fc-2 



V3 + 2w 

If fc 7^ 2, we find that in the Einstein frame, when the Universe isotropises 

oc t 2 '^ 



e~ il oc to^) 1 

with k > 2 so that the scalar field is real. It follows from go) that t oc ^-2)/(fc-i) anc j 
near the isotropic state, the scalar field and the metric functions behave in the Brans- 
Dicke frame like : 

OC ^/C 1 -*) 
g-H K ^l + fc+2w )/((fc-2)(fc-l)) 



When fe = 1, one has i oc e* and thus oc e 2t and the Universe in the Brans-Dicke 
frame tends to a De Sitter model. 

When k — 2, I is zero and the scalar field asymptotically tends to a constant. Con- 
sequently, t —> i and the Universe tends to a De Sitter model in the Einstein and 
Brans-Dicke frames. 

These results agree with those of |12| and are illustrated on figures |2] and [5] for the 
Bianchi type / model. The first one is such as I 2 < 3 and shows that the metric func- 
tions derivatives converge to a same behaviour whereas the metric functions grow at 
the same rate. The second one is such as I 2 > 3 and the Universe does not isotropise : 
the metric functions derivatives never converge to a common behaviour whatever t and 
the metric functions deviate from each other. Some similar figures may be obtained for 
the three others scalar-tensor theories that we consider in this section. 





10 20 30 40 50 60 70 



Fig. 2 — Logarithm of the metric functions (first graph) and their derivatives (second graph) with respect to the proper 
time and luq — 2.3 and k — — 1. Then £ 2 — 1.18 and the Universe isotropises : the metric functions derivatives tend to 
a common behaviour,the metric functions growing at the same rate. 



- Brans-Dicke theory with an exponential potential 

This theory is defined by G _1 = (f>, Q — o>o and U = e k< ^ and is thus such as U = 
e fc00-2 i j s men a function of the scalar field : 

k(f>- 2 



t = 



V3 + 2cj 
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Fig. 3 — Logarithm of the metric functions (first graph) and their derivatives (second graph) with respect to the proper 
time and lDq — 2.3 and k — — 3. Then £ 2 — 3.28 and the Universe does not isotropise whatever the length of the time 
integration : the metric functions derivatives do not tend to a common behaviour and the metric functions deviate from each 
other. 



In the Einstein frame, we calculate in the neighbourhood of the isotropic state that the 
scalar field behaves like : 

</) = 2/(fc-0 o e 4n /( 3+2 ^) 

where <f>o is an integration constant. It thus tends exponentially to the constant 2/k 
when fl — > — oo. Consequently £ tends to zero, t — ► t and the Universe tends to a De 
Sitter model in the Einstein and Brans-Dicke frames. 

- String theory with an exponential potential 

This theory is defined by G _1 = e~^, Hi = <pe^^ and U — e k< ^ . Consequently, U — 

e (k+2)<f> an( j j t comes 

t _ (2 + fc)0 

In the Einstein frame and in the vicinity of the isotropy, we deduce from Q8t that 

ExpIntegralEi(-4>)) - 3/0 = 2(2 + k)Q + (j) 

where ExpIntegralEi(—(f)) — e~ x x~ 1 dx. We deduce that if k > —2 (respec- 
tively k < —2), 51 — > — oo when (/> — > 0+ (respectively — > 0_). Then a power series 
of this equation about <fr = 0, learns us that 



2(2 + k)n + 4> 

Consequently £ — > like —3 [2(2 + fc)f2]~ , and the Universe tends to a De Sitter 
model in the Einstein frame. We also find that f oc i and thus in the Brans-Dicke frame, 
<fi oc t^ 1 oc t^ 1 and the Universe again tends to a De Sitter model. It is the same for 
the special value k = — 2. 

- String theory with a power potential 

This theory is defined by = e"^, uj = and U = (j> k and thus U = e 2 ^(j) k . It 
comes that 

g _ fc + 20 
v/3 + 20e~* 
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and we asymptotically find that in the Einstein frame the scalar field is solution of 



e k/2 ExpIntegralEi(-- -</>) + - In - n , = 20, + <p 
2 k k + 2<p 

Consequently 57 — > — oo when </> with k > and a power series of this equation 
about = shows us that 

It implies that near the isotropic state £ —> k/y/3 and in the Einstein frame 

0->f- 2fc_1 
e -n ^ t 3*- a 

Consequently t ~ i and it comes in the Brans-Dicke frame 

These results show the importance of the power or exponential law of the proper time 
as an asymptotic form for the metric functions in the vicinity of the isotropy in both 
Einstein and Brans-Dicke frames. In the following section, assuming these forms of 
funtions for e , we examine the asymptotical behaviours of the gravitational function 
G and potential U of the scalar field when the Universe isotropises. 

4.2 Isotropisation of a scalar-tensor theory when e n is a known 
function of the proper time 

The above theories show that when the Universe becomes isotropic, the metric 
functions of the Brans-Dicke frame asymptotically behave like some power or expo- 
nential functions of the proper time, typical of the solutions found without a scalar field 
or with a cosmological constant. In fact they are among the most usually widespread 
asymptotic behaviours in the literature (De Sitter model, FLRW models with a perfect 
fluid, etc). Reversing the steps of the previous section, we are going to determine the 
asymptotic forms of the potential and the gravitation function near the isotropic state 
when the metric functions of the Brans-Dicke frame tend to tr or e kt . We will thus 
get the conditions allowing to G and U to be bounded when the Universe isotropises, 
which currently seems to be the case, and will compare them with those allowing the 
isotropisation of the Universe. For this purpose, we will use the following algorithm : 

1 . For each form of f2(t), we will suppose that £ tends to a vanishing or non vanish- 
ing constant what implies respectively that e~ n tends to t or to an exponential 
of t. 

2. We then calculate G^ 1 (t, t) using (EH . 

3. We deduce from it i(t) using (gOJ and thus G _1 (f) and U(t) using J42l 
Employing this method, we obtain the following results : 

- When e~ n — > i k and the Universe isotropises such as £ — > cte ^ 0, G _1 — > 

p(l-fc£ 2 )/(£ 2 -l) an( j fj = t-2[(.l+k)l 2 -2]/(l 2 -l)^ 
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- When e~° — > i k and the Universe isotropises such as I — > 0, G" 1 — > i~ 2 and 
[7 — > t . It is in accordance with the previous point. 

- When e — > e fe * and the Universe isotropises such as t — > c<e 7^ 0, G -1 — > 

e 2fcf 2 /(i-« 2 )* an d [7 -> e 2fc« 2 /(i-£ 2 )f 

- When e _n — * e fc * and the Universe isotropises such as I — > 0, G -1 — * G^ 1 and 
U —> Uq, Gq 1 and Ob being two constants. It is in accordance with the previous 
point. 

These results are in agreement with those of the preceding section and generally stable 
as long as e~ n and e approach their asymptotic behaviours (in power or exponen- 
tial of the proper time) more quickly than respectively t^ 1 and t _1 . We discuss their 
significance in the following section. 



5 Discussion 

In this work, we explained in section|2]the main steps to obtain the ADM Hamil- 
tonian of the General Relativity with a minimally coupled and massive scalar field. We 
then showed in section |3]how to use this formalism to find, for the homogeneous but 
anisotropic flat Bianchi type / model, the dynamics of the metric functions according 
to the initial and final values of the Hamiltonian as well as the sign of the potential. 
We now comment these results. In what follows we will call "Big-Bang singularity" 
a pointlike singularity characterised by the vanishing of the three metric functions. A 
"pancake singularity" or "cigar singularity" will apply, respectively to the cases where 
one or two metric functions vanish. 

From the table El we obtain the following results when the sign of the potential does 
not vary : 

- A metric function can have a maximum (minimum) only in presence of a nega- 
tive (positive) potential. Moreover, all the conjugate momenta p$j can not have 
the same sign and then the condition C\ can not be true for all the metric func- 
tions : it follows that at least one metric function will not have any extremum. 

- Since C\ can not be checked for all the metric functions, when the Hamiltonian 
is positive, the three metric functions can be increasing together at late times, but 
not decreasing. All types of singularity, Big-Bang type, pancake type or cigar 
type are possible at early time. 

- Since C\ can not be checked for all the metric functions, when the Hamiltonian 
is negative, the three metric functions can be decreasing together at late time but 
not increasing. The singularity if it exists will only be of pancake or cigar type 
at early time. 

- From Collins and Hawking isotropy definitionl 13 1, the Universe will isotropise 
when 57 — ► —00. It corresponds to the vanishing of df3±/dt oc e 3n and is neces- 
sary to the convergence of (dgij /dt) / gij to a common value allowing to observe 
the same Hubble function on all the directions. We thus easily deduce, that Uni- 
verse isotropisation will arise at late (early) t times if the Hamiltonian is positive 
(negative). 

As written in the first point, as long as the potential has a constant sign, the metric 
function can have one and only one extremum. This is also the case when we consider 
a flat or open FLRW model with trace-free matter, (f> finite and co^ > as shown in 
II 1 41 - In this paper it is also proved that flat FLRW models can only contain a single 
minimum whereas here, we shown that a single maximum is also allowed for a nega- 
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tive potential. When the sign of the potential varies, the results of table |2] are always 
true but Him and H en d define the intervals of values of the Hamiltonian for which the 
sign of the potential is constant. 

Using a conformal transformation of the metric whose rules are given by (1401421 and 
the results about the Universe isotropisation obtained within the framework of a min- 
imally coupled and massive scalar field, we studied in section |4] the isotropisation of 
several scalar-tensor theories, spread in the literature, for the Bianchi class A models. 
We obtained the following results in the Brans-Dicke frame : 

- Brans-Dicke theory : 

- With a power potential of the scalar field U = <jr, a necessary condition for 
the isotropisation will be that (k — 2) 2 (3 + 2tI>o)~ 1 < 3 or < 1 according 
to the absence or the presence of curvature. So that the scalar field is real, 
it is necessary that k > 2. Then when k > 2, the metric functions tend to 
£i+k+20o)/((fc-2)(k-l)) and the Universe is expanding if I + k + 2us Q > 0. 
Moreover, this expansion will be accelerated if (k — 2) 2 (3 + 2ujq)~ 1 < 1. 
Hence in presence of curvature, the Universe expansion is always accelerated 
when it isotropises. The potential U tends to vanish as t 1 ^ , thus being able 
to explain the smallness of the cosmological constant at late times. However, 

_ 2 

the gravitation function G diverges as t k - x in disagreement with its small 
observational value. If k — 1, the Universe isotropises without condition to 
a de Sitter model, the gravitation function tends to vanish but the potential 
diverges. If k = 2, the Universe isotropises without condition to a de Sitter 
model and U and G tend to some non vanishing constants. 

- With an exponential potential of the scalar field U — e k< ^, the Universe isotropises 
without condition to a De Sitter model and is expanding. U and G tend to some 
non vanishing constants. 

- String theory at low energy : 

- With a power potential of the scalar field U = cj) k , a condition necessary 
to the Universe isotropisation will be that k 2 < 9 or < 3 according to the 
absence or the presence of curvature and k > 0. The metric functions tend 
then to P k and the expansion of the Universe is thus accelerated ifk 2 < 3, 
which is always the case in presence of curvature. G tends to a non vanishing 
constant and U vanishes like iT 2 , explaining the smallness of the cosmological 
constant at late times. 

- With an exponential potential U = e k< ^, a necessary condition for isotropisa- 
tion will be that k > —2(k < —2) if the scalar field is positive (negative). 
The metric functions tend then to a De Sitter model and U and G to some non 
vanishing constants. 

These applications clearly show, contrary to what occurs for a minimally coupled 
scalar field, that the Universe isotropisation does not lead systematically to an expand- 
ing Universe. This comes owing to the fact that the conformal transformation does not 
preserve the signs of the metric functions derivatives : an expanding Universe in the 
Einstein frame may be contracting in the Brans-Dicke frame. If these two theories are 
among the most studied in the literature, power or exponential laws of the proper time 
for the metric functions are also the most widespread. This is why it appeared signifi- 
cant to us to study them. Hence, the De Sitter model is often used to solve the flatness 
problem and lots of FLRW models solutions with a dust or a radiative fluid gener- 
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ally behave like powers of the proper time. In 1 15 1 where the asymptotic dynamics of 
the scalar-tensor theories with a perfect fluid are studied for the FLRW models, the 
late time solutions of the metric functions when they converge to General Relativity, 
give place to power or exponential laws of the proper time. It is also the case in 1161 
where one looks for scaling attractors likely to produce an accelerated expansion. It 
thus strongly justifies the study of these forms of metric functions as the outcome of 
isotropisation process. One obtains the following conditions in the Brans-Dicke frame 
ensuring the non divergence of the potential and the gravitation function when the Uni- 
verse isotropises and e~ n tends to t k or e kt : 

When the Universe isotropises, is expanding and that asymptotically the metric func- 
tions behave like a power of the proper time i (k > 0), it is necessary and sufficient 
so that the potential and the gravitation function do not diverge asymptotically that 

- when t 2 < 1, ^ < £ 2 < y-2-£ ; thus k > 1 and the Universe expansion is always 
accelerated. 

- when i 2 > 1, i.e. only for a flat Bianchi type I model since isotropisation does 
not occurs for a model with curvature in such a case, ttt < I 2 < \ '■ thus k < 1 
and the Universe expansion is thus decelerated 

If the metric functions behave asymptotically like an exponential of the proper time e kt 
(k > 0), the potential and the gravitation function do not diverge only iff. — » 0. 

This important result shows that isotropisation andfiniteness of the potential and grav- 
itation function only occur when the Universe tends to a De Sitter model with I — > or 
when the metric functions tend to a power law of the proper time representing an ac- 
celerated Universe with a (vanishing or not) curvature or a decelerated flat Universe. 
Hence, the expansion being accelerated today, this result is clearly in accordance with 
the presence of at least a small curvature in our present Universe. 
One can also obtain some results about the No-Hair theorem which states that in pres- 
ence of a cosmological constant, the Universe tends to a De Sitter model. With regard 
to General Relativity with a massive scalar field, we showed in |17| that this result 
could be generalised to any theories such as I — > (when the assumption 4 of the 
section|4]is respected) which thus include the particular case U — constant for which 
£ = strictly. For the scalar-tensor theories, the fact that I tends to zero is not sufficient 
any more to imply the convergence of the Universe to a De Sitter model because the 
asymptotic behaviour of the metric functions also depends on that of <f>. On one hand, 
if the potential U is a constant, I is exactly zero and the scalar field tends to a con- 
stant, implying that the Universe behaves like a De Sitter model. On the other hand, if 
the potential U is a constant but that I tends to a non vanishing constant, G _1 — > t 2 
and e _sl — > t^ 2 ^ 1 . Consequently, the presence of a cosmological constant in the 
Brans-Dicke frame does not imply necessarily the convergence of the Universe to a De 
Sitter model contrary to what occurs in the Einstein frame . 

All these results were established thanks to the ADM Hamiltonian formalism. It is 
a powerful tool but which unfortunately remains difficult to use for the Bianchi class 
B models. It would thus be interesting to deal with this class of models by using the 
results of 1 7 1 in which the divergence theorem is redefined. 
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